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ABSTRACT 

A technique involving 55 Fe X-rays provides a straightforward method to measure the 
response of a detector. The detector's response can lead directly to a calculation of the 
conversion gain (e _ ADU -1 ), as well as aid detector design and performance studies. We 
calibrate the 55 Fe X-ray energy response and pair production energy of HgCdTe using 8 
HST WFC3 1.7 /im flight grade detectors. The results show that each Ka X-ray generates 
2273 ± 137 electrons, which corresponds to a pair-production energy of 2.61 ± 0.16 eV. The 
uncertainties are dominated by our knowledge of the conversion gain. In future studies, we 
plan to eliminate this uncertainty by directly measuring conversion gain at very low light 
levels. 

Subject headings: Astronomical Instrumentation, Infrared Detectors, Conversion Gain, In- 
terpixel Capacitance 
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1. INTRODUCTION 

Conversion gain, g c (e _ ADU -1 ), is a fundamen- 
tal parameter in detector characterization that is 
used to measure many detector properties, includ- 
ing read noise, dark current, and quantum efficiency 
(QE). Measuring these parameters to a higher preci- 
sion is becoming more important as the demand for 
low-signal observations increases and the scientific 
requirements evolve in complexity. 

The upcoming Joint Dark Energy Mission 
(JDEM), for example, requires the precise measure 
of the dark energy equation equation of state, w 



and 



This measurement is derived from the 



luminosity distance versus redshift relationship of 
Type la supernovae out to high redshifts. To ob- 
tain the necessary accuracy, the mission requires de- 
tector characterization measurements to within 2% 
(jSchubnell et alJll^ . 

A traditional method for measuring conver- 
sion gain is the photon transfer technique 
(jJanesick et al.1 [19871 ) . This method, however, is 
limited in that it does not account for charge cou- 
pling between pixels. Charge coupling attenuates 
the photon shot noise of individual pixels, resulting 
in an underestimated variance and an overestimated 
gain. Figure [1] illustrates the coupling of the central 
pixel with the four nearest neighbors, both for (a) 
symmetric and (b) partially-symmetric cases. The 
degree of coupling is given by coupling coefficients, 
a and (3. 
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Interpixel capacitance (IPC), one of the more 
dominant charge coupling mechanisms, is caused 
by small amounts of stray capacitance. Early 
simulations first pre dicted the effects of IPC 
(Kavadia s et al.lll994l ). Later studies measured the 
degree of couplin g using various m ethods, including 
cosmic ray hits (iFiger et alj|2004|). e xternally cali- 
bra ted capacitors (iFinger et al.ll2006f). an d hot pix- 
els (|Bai et al.ll2007t ). iMoore et all (|2006l) provide a 
detailed mathematical discussion. These works all 
indicate that ignoring IPC effect within infrared de- 
tectors ca n result in a n over estimate of gain by up 
to 20%. iDorn et al.l (|2006j ) show these effects to 
be as a high as a factor of 2 within the fully de- 
pleted silicon-based hybrid CMOS arrays, such as 
the Rockwell HyViSI H2RG arrays. 

While the conversion gain is ultimately measured 
via standard stars once the in strument has bee n 
launched and commissioned (e.g. lReach et al.l2 005). 
detectors must meet specifications and science mis- 
sions must be planned much earlier. An overesti- 
mate of the gain can significantly distort the assess- 
ment of the detector and mission potential. IPC 
becomes more noticeable in near-infrared detectors 
as pixel size decreases. Better methods of measur- 
ing g c therefore have broad applicability to astro- 
nomical detector characterization. For example, one 
popular method relates the strength of the inter- 
pixel coupling and, subsequently, the "true" vari- 
ance to the autoco rrelation function of the detector 
(jMoore et al. 2006). We discuss this method further 
in Section [3~2l 

A technique involving 55 Fe X-rays provides an al- 
ternative, straightforward method to measure the 
conversion gain. This technique quantifies the nec- 
essary energy to liberate a single electron-hole pair 
in the semiconductor, known as the pair creation 
energy. If the X-ray energy is known, the number 
of electrons liberated by each X-ray photon follows. 
The conversion gain is then derived from the ob- 



served instrumental counts (ADU) produced by a 
single X-ray. Already, this is a popular technique 
within the CCD community, where each 55 Fe Ka 
X-ray pho ton is known to li berate ~1620 electrons 
in silicon (jFraser et al.lll994t ). 

This absolute calibration offers additional bene- 
fits that aid detector design and performance stud- 
ies. For example, the calibrated response provides 
a fundamental standard for measuring drifts and/or 
changes in a detector's system gain over time. Addi- 
tionally, if the conversion gai n is given by a reliab le 
alternative method, such as iMoore et al.l |2006)'s 
autocorrelation technique, the 55 Fe X-ray response 
serves as a direct probe of the internal quantum effi- 
ciency (QE), also known as the charge collection effi- 
ciency In other words, this technique measures the 
detector response to individual X-rays so that the 
QE does not include any contribution from external 
factors such as the anti-reflective coating or long- 
wavelength photons that are not absorbed. While 
not commonly used by the astronomical community, 
the internal QE is an important measurement to de- 
tector manufacturers in their quest to build more 
efficient detectors. 

A semi-empirical trend exists between the 
pair creation energy and semiconductor bandgap 
(|Alig fc Bloo"nl[T975h : 

e = 2.73^ g + 0.55, (1) 

where e is the pair-production energy (in eV) 
and Eg is the band-gap energy (in eV). An ex- 
act pair-production energy, however, is not known 
for the commonly used mercury-cadmium-telluride 
(HgCdTe) infrared detector arrays. The limit- 
ing factor has been the cadmium-zinc-tellruidc 
(CdZnTe) substrate found on many HgCdTe de- 
tectors, which absorbs the 55 Fe X-rays before they 
reach the depletion region of the detector. The in- 
frared detectors for the HST Wide-Field Camera 3 
(WFC3) are not limited in this respect because they 
are all substrate-removed. 

In this paper, we provide calibrated measure- 
ments of the 55 Fe X-ray response of 1.7 /jm HgCdTe 
detectors. Section [2] presents our 55 Fe data. To con- 
vert the X-ray intensity from ADU to electrons, we 
establish the "true" variance and gain of a pixel 
using classical propagation of errors. Section [3] dis- 
cusses the mathematical framework for this method. 
Monte Carlo simulations validate the fidelity of this 
technique. Section [4] presents the results. We 
find that in response to the 55 Fe X-rays, these 
1.7 ^m HgCdTe detectors yield 2273 ± 137 elec- 
trons, which corresponds to a pair-production en- 
ergy 2.61 ± 0.16 eV. We describe in detail the vari- 
ous systematic errors inherent within this measure- 
ment. 

2. TEST METHODOLOGY AND DATA 

We have collected 55 Fe X-ray data from eight 
Teledyne Imaging System's (TIS) 1.7 H1R 
flight grade WFC3 detectors, which are molecular 
beam epitaxy (MBE) grown and CdZnTe substrate- 
removed. These detectors have a read noise of —15 
e~ rms per read (—21 CDS) and a dark current 



of —0.02 e~ s~ 1 at 145 K. All data were collected 
at the NASA Goddard Space Flight Center (GSFC) 
Detector Characterization Laboratory (DCL). The 
55 Fe X-ray source was placed at a normal incidence 
1.5 cm from the center of the 1.8x1.8 cm FPA. The 
data were all collected in 2008 and have since been 
archived. 

For each detector, we obtained 10 exposures, each 
consisting of 6 correlated double samples (CDS) 
taken up-the-ramp, where each CDS pair is com- 
posed of the th and n th frames up the ramp, where 
n is a number between 1 and 6. We performed a 
standard dark subtraction with median frames cre- 
ated from the 10 exposures. We include plots for 
only FPA 152, which is a typical representation of 
the other detectors we tested. We discuss any dif- 
ferences observed in the other detectors. 



Figure 2(a) shows a histogram of individual pixel 
intensities (in ADU) from the 10 exposures of FPA 
152. Key features include the dominant Ka peak, 
the less dominant K/3 peak, and a low-intensity 
plateau (that extends down to zero signal). The 
two peaks show broadening likely due to a combi- 
nation of several effects. Off-centered X-ray hits and 
charge diffusion, which occurs on the 0.1-pixel level, 
cause broadening to the left. A changing band-gap 
energy with penetration depth causes broadening 
to the right. Read noise, pixel-to-pixel gain varia- 
tions, and the intrinsic fano factor 8 cause overall 
broadening. We fit the peak of each profile us- 
ing IDL's GAUSSIAN procedure, which performs 
a combined gaussian and polynomial fit. The low- 
intensity plateau results mostly from charge diffu- 
sion and non-centered hits. 

The single pixel histogram in Figure 2(a)| how- 
ever, ignores charge coupling effects that cause some 
of the signal from the X-ray to appear in neighbor- 
ing pixels, such as in Figure [TJ The histogram in 
Figure [3] shows the fraction of the total charge in 
FPA 152 that couples to neighboring pixels (where 
the total charge is defined as the intensity in the 
central pixel and four nearest neighbors). The plot 
reveals a 1-2% coupling in all primary directions. 
While there is a distinct symmetry in both the hor- 
izontal and vertical directions, the horizontal cou- 
pling tends to be larger than the vertical coupling. 
This partial symmetry suggests a non-symmetric 
multiplexer associated with the horizontal readout 
mechanism. Overall, all nearest neighbors (solid 
line) show an average coupling of —1.5% while a 
distribution of random pixels (dotted line) shows, 
as expected, no coupling to the X-rays. 

Adding the nearest neighbor flux presents the 
most straightforward method by which to restore 
the lost charge. We do not perform aperture pho- 
tometry. While coupling is observed in other pixels, 
such as diagonal and second-neighbor pixels, Monte 
Carlo simulations reveal this effect to be negligible 

8 The Fano Factor adjusts for additional line broadening 
due to the statistical nature of the loss of energy to thermal 
excitation in the lattice. An accurate measure of the Fano 
factor requires detailed modeling of proprietary data from 
TIS. 



on the 1% level. The additional read noise from 
these pixels is greater than the actual signal. We 
conclude that the additional noise from these pix- 
els would be more detrimental to the total signal 
than the small amount of charge lost when adding 
only nearest neighbors. We consider only hits that 
are isolated (defined within a 5 x 5-pixel area) and 
well-centered (defined as bei ng associated with the 
peak intensity in Figure 2(a) as opposed to the low- 
intensity plateau). This thresh old is indicated by 
the vertical line in Figure 2(a) In this case, "well- 
centered" is a somewhat arbitrary definition, and 
we find that shifting the threshold can alter the re- 
sults by tens of ADU. We discuss the implications 
of this effect in our discussion of systematic errors 
in Section P 



Figure 2(b) plots a histogram of only the selected 



pixel intensities with the nearest neighbor charge 
restored. Expectedly, the peak of this distribution 
shifts to a higher intensity, consistent with the frac- 
tional coupling observed in Figure[3J The additional 
read noise from the neighboring pixels blends the 
Ka and K/3 peaks. We again take the peak by fit- 
ting the distribution with IDL's GAUSSIAN rou- 
tine. The peak in this case represents the detec- 
tor's response (in ADU) to a single 55 Fe X-ray, for 
which the energy now includes a weighted contri- 
bution fro m the K(3 line. Given a Ka-to-K/3 peak 
ratio of 7:l (|Scofieldlll974f ). we approximate the X- 
ray energy as a weighted average: 5.9375 keV. 

3. CALCULATING GAIN VIA PROPAGATION OF 
ERRORS 

The pair-production energy of the detectors ul- 
timately requires the total intensity of the X-rays 
to be in electrons. The conversion of the detector 



response in Figure 2(b) from ADU into electrons, 
therefore, must be calculated with an independent 
measure of the conversion gain. The photon trans- 
fer method is insufficient because charge coupling 
attenuates the photon shot noise and decreases the 
measured variance. Here, we explore a technique to 
establish the "true" variance via the classical prop- 
agation of errors. The technique has the advantage 
of allowing for a direct measurement of the covari- 
ance matrix from the data, which can serve as a 
diagnostic of other correlations that may be present 
in addition to IPC. We find this method works well 
at all signal levels an d is consistent with the results 
of lMoore et al.l (|2006| ). We therefore implement this 
method throughout our analysis. 

3.1. Mathematical Preliminaries 

We begin by considering the simplest possible ex- 
ample of small symmetric crosstalk to the four near- 
est neighbors. (A partially symmetric crosstalk con- 
sistent with Figure [3] is more complicated and con- 
sidered in Section [Aj] ) As we discuss in Section^ 
we ignore diagonal and second-neighbor pixels be- 
cause coupling effects are negligible on the 1% level 
for our detectors. For larger IPC, other neighbors 
must be considered. 

The detector is uniformly illuminated and only 
Poisson noise is present (i.e. no read noise). This 



model is true for high signal levels for which shot 
noise in the signal dominates. Crosstalk is assumed 
to be symmetric and to only affect the nearest neigh- 
bors. A 3 x 3 pixel subarray is considered: 



s si s 2 

S3 S4 S 5 

s 6 s 7 s$ 



(2) 



The interpixel charge coupling can be represented 
by a detector point spread function (PSF) that, 
when convolved with the subarray, "blurs" the 
charge between pixels. This kernel is represented 
as: 

/ a \ 
PSF = [a 1 - 4a a , (3) 
V a / 

where a is the coupling coefficient. In other words, 
a is the fractional crosstalk between the central 
pixel and its nearest neighbors. The "true" signal, 
s, in a pixel is therefore a linear combination of the 
measured signal in the central pixel of interest, s 4 , 
and its neighbors, 



i=0 



(4) 



where ai is some fraction of the signal in pixel Sj. 
In fact, Equation 0] represents the convolution of a 
kernel that undoes the blurring caused by the PSF. 
In other words, a must be the inverse PSF, 







PSF * PSF = 1 , 

V J 



(5) 



where * is the convolution operator. For the PSF 
of Equation [3l the inverse PSF is, 



2a -l + 2a 2a 

-l + 2a -5 + 1 + 2a -1 + 2a 



PSF 



1 - 9a + 18a 2 



2a 



-l + 2a 2a 
(6) 

By the associative property of convolution, 



f*(g*h) = (f*g)* h, 
Equation [¥] can be written as, 
1 



(7) 



s 4 



(s 4 + a (— si — s 3 - 5s 4 — S5 — s 7 



l + 9a(-l + 2a) 
+2a (s + si + s 2 + s 3 + s 4 + s 5 + s@ + s 7 + s 8 ))) 

The "true" variance of the pixel can be shown 
using the classical propagation of errors, 



(8) 



8 8 



\ -\ - ds 4 ds 4 



j=0 i=0 



dsi dsj 



(9) 



where Cij is the covariance between pixels i and j. 
The covariance matrix is given as 

1 N 

c v = j^jY,^- j ^~j)> ( 10 ) 



where N is the number of pixels. iMoore et al.l 
(2006) show that to bring a one percent correlation 
signal up to the noise level, the number of samples 
must satisfy the equation 

iV>4j. (11) 
or 

A coupling coefficient of ~1% requires > 10 6 sam- 
ples. 

The covariance matrix can be found experimen- 
tally (see Section [A.2[) with > 10 6 independent in- 
tegrations, but our archival data set is limited to 
< 10 integrations for each FPA. To obtain the neces- 
sary number of samples in a realistic period of time 
requires a unique detector readout mode capable of 
sampling only small windows of pixels. While such 
readout modes are currently being implemented for 
projects such as the James Webb Space Telescope 
(JWST), these readout modes are unavailable for 
HST because they are not required by the planned 
science. 

Alternatively, the covariance matrix may be de- 
rived from theory under the assumption that noise 
in neighboring pixels is correlated solely with shot 
noise in the pixel of interest. The theoretical covari- 
ance matrix is given as: 
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Equation[T2]requires knowledge of only the coupling 
coefficient, a, and the measured noise, a. We use 
Equation 1 121 throughout our analysis. 

Substituting Equations [8] and [12] into Equation [9] 
gives the "true" variance in the pixel of interest, 

a 2 S4 « (1 + 8a + 52a 2 + ...)<, (13) 

where <r 2 4 is the measured variance. IMoore et al.l 
(2006) derived the same results to first order in a 
using Fourier methods. 9 

3.2. Validation 

Monte Carlo simulations validate the propagation 
of errors technique given by Equation [T^l Each 
simulation models charge coupling in a uniformly 
illuminated 1024 x 1024 detector array. Illumi- 
nation levels range in brightness from 10 < s < 
10 5 e~ pixel -1 , where s is the mean integrated flux 
per pixel. The shot noise is distributed with a Pois- 
son distribution. Charge coupling is simulated via 
convolution of a PSF given by Equation [3] with a 
coupling coefficient a — 1.5%. A Gauss-normal dis- 
tribution of read noise, <7 rea d = 15 e~, is added last. 
We choose these parameters to be consistent with 
the actual data observed in Section [U 

9 They did not explicitly write out higher order terms. 



The "true" variance of the data when no charge 
coupling effects are present is given as: 

cfrue = S + ofead ( 14 ) 

The simulated charge coupling attenuates the vari- 
ance. We implement three different techniques to 
recover the "true" variance: (1) a l+8a approxi- 
mation, (2) the au tocorrelation technique given by 
IMoore et alj (|2006), and (3) the propagation of er- 
rors technique described above in Section [3. II 

The first method attempts to restore the variance 
by multiplying by a 1 + 8a coefficient. This co- 
efficient accounts for the ~ 1 — 8a variance loss 
gi ven by a small a a pproximation in equation 28 
of IMoore et al.l (|2006D . The technique is commonly 
used within the near-infrared community because of 
its straightforward implementation. The coupling 
coefficient, a, is typically measured during a dark 
exposure by utilizing the single pixel reset feature 
of the H2RG multiplexer to set a si ngle pixel to a 
volta ge different from its neighbors (|Seshadri et al.l 
2008). The crosstalk induced by IPC on neighbor- 
ing pixels is observed directly. 

The auto c orrela tion procedure described by 
IMoore et al.l (|2006l ) relates the "true" variance to 
the autocorrelation matrix: 

4 ue = 2^(Z> 2 [*>j] + 2 E^<jW + m+ 

2j2D[iJ]D[i,j+l}), (15) 

where D is the difference of two equally illumi- 
nated flat fields and N is the number of pixels. 
Equation [15] makes a small a approximation and 
considers all contributions from pixels beyond the 
four nearest neighb ors to be negligible. In addition, 
IMoore et alj (]2006f ) assume IPC to be the dominant 
coupling mechanism. Other coupling mechanisms, 
however, typically are present, such as correlations 
in the the photon arrival (e.g. Bose=Einstcin or 
Hanbury-Brown-Twiss) and carrier diffusion (e.g. 
carrier-carrier interaction). Despite these assump- 
tions, this variance estimator i s typically favored 
by the near-infrare d community (|Brown "etaL|[200l 
iFinger et alll200ll . 

Finally, method three implements the propaga- 
tion of errors technique described in Section 13.11 
The "true" variance is given by Equation [T3l The 
covariance matrix in Equation [12] directly measures 
contributions from all additional coupling mecha- 
nisms that may be present. Additional data could 
potentially resolve the different coupli ng mecha- 
nism contributions. Nonetheless, like IMoore et al.l 
(2006), we assume IPC to be the only coupling 
mechanism for the purposes of this paper. We sim- 
ulate our real data by creating only 10 integrations. 
Since the propagation of errors technique requires 
10 6 independent integrations to calculate the covari- 
ance matrix experimentally, we instead use the the- 
oretical form given by Equation 1 12[ where the cou- 
pling coefficient, a, and the noise, a, are measured 
directly. Section 13.31 describes the measurement of 
the coupling coefficient in more detail below. 



Figure 0] plots the ratio of the measured to the 
true variance versus signal for each of the above 
methods. As we predict, no correction results in an 
underestimated variance. The rise at small signals 
is due to the dominating read noise component. At 
high signal levels (> 500 e _ ), the 1 + 8a approxi- 
mation effectively estimates the variance, although 
there is a disagreement of ~1% due to the small 
a approximation. At low signal levels, however, 
this technique is ineffective because the dominant- 
ing read noise does not couple between pixels. The 
autocorrelation technique restores the "true" vari- 
ance to within 1% at all signal levels. This method, 
however, does not fail at low signal levels because 
read noise does not contribute to the autocorrela- 
tion terms. 

The propagation of errors technique is most effec- 
tive at restoring the "true" variance at all signal lev- 
els. Again, there is no read noise contribution to the 
autocorrelation coefficients. This technique is not 
as straightforward as the autocorrelation technique 
because it requires a direct measurement of the cou- 
pling coefficient, but it allows us to restore the vari- 
ance most effectively. For this reason, we implement 
the propagation of errors technique throughout the 
rest of this paper. 

3.3. Measuring the Coupling Coefficient, a 

The propagation of errors technique (Equa- 
tion I13p requires knowledge of the coupling coef- 
ficient, a. Unlike the simulations in Section 13.21 
where the coefficient is a known input, the actual 
coupling c oefficient must be m easured directly from 
the data. I Mo ore et alJ (|2006l ) show that a can be 
extracted directly from the autocorrelation matrix 
of uniformly illuminated exposures. The ratio, r, of 
the nearest neighbor and unshifted autocorrelation 
coefficients is given as: 



2a(l - 4a) 
(f - 4a) 2 +4a 2 ' 



(16) 



Table [T] lists the autocorrelation matrix of FPA 
152. The other detectors yield similar matrices. 
The autocorrelation matrix is computed from the 
differenced image of two identically illuminated ex- 
posures. To account for l//-noise variations, we 
averaged differenced images computed from pairs 
of flat-field exposures with varying delay times be- 
tween exposures. We assume symmetry in both the 
horizontal and vertical directions so that the auto- 
correlation of each pixel with itself (i.e. unshifted) 
is given by the upper left ([0,0]) matrix value. 

All values are normalized to the unshifted auto- 
correlation coefficient. The nearest neighbors ([0,1] 
and [1,0]) have correlations that are typically be- 
tween 1-5%. Even diagonal and second neighbor 
pixels show some correlation although the degree is 
much less than the nearest neighbors. Solving for 
a, Equation [16] indicates the coupling coefficients 
are typically between 1-2%. Figure [5] plots the cor- 
rected conversion gain versus signal and compares 
it to the uncorrected values for FPA 152. As pre- 
dicted by the Monte Carlo simulations, the uncor- 



rected gain is typically overestimated by ~ 13% for 
coupling coefficients of ~ 1.5%. 

3.4. Measuring the Conversion Gain, g c 

The relevant conversion gain, g c , should corre- 
spond to signal levels equivalent to the 55 Fe X-ray 
intensity, which is highlighted by the vertical line in 
Figure [5] (~800 ADU). The figure reveals, however, 
that the signal levels of the gain data are typically 
higher. While an obvious solution would be to re- 
take the necessary data, our analysis is limited to 
archival data. Instead, we attempt to estimate the 
conversion gain at low signal levels and minimize 
the associated systematic error. 

Figure [5] shows a correlation between the con- 
version gain and signal level. This trend arises 
from the relationship between the capacitance of 
the diode and bias voltage , which is well understood 
for HAWAII devices (e.g. iFinger et al.ll2006fl . The 
capacitance increases with signal level because the 
voltage, and hence the depletion width of the pn 
junction, across the diode decreases. For low signal 
levels (i.e. small capacitances) an electron generates 
a larger voltage swing and the gain (in e~/ADU) is 
smaller. 

To better understand the gain behavior of these 
devices at low signal levels, we performed alter- 
nate gain measurement techniques, including an 
analysis of the photon transfer of warm pixels 
in dark ramps and sol ving for the gain in dark 
ramps via equation 1 of iRauscher et alJ (2007). We 
also compared our results to ultra-low light level 
data collected for similar FPAs undergoing testing 
for the JWST Near-Infrared Spectrograph project 
(Rauscher, priv. comm.). These results all confirm 
the gain decreases with signal level and that the 
decrease is smooth and continuous. 

We therefore choose to fit a curve to the high sig- 
nal gain data and extrapolate to lower signal levels. 
The choice of curve to fit is unclear because the true 
behavior of the gain is not well described at these 
low signal levels , where the data tend to be nois- 
ier. Figure 5(a) s uggest s that an exponential fit is 
sufficient. Figure 5(b) however, indicates that a 
polynomial provides a better fit. We find the likely 
gain is bound by the extremes of these two curves 
and take the conversion gain to be an average of 
the two. Table [2] lists the "true" gain for all the de- 
tectors. We discuss the systematic error associated 
with this method in Section [ 



4. RESULTS 

4.1. Pair- Production Energy 

We calculate the pair-production energy by divid- 
ing the X-ray response (in e~) into the energy of a 
single 55 Fe X-ray, which includes a weighted con- 
tribution from the blended K/3 peak. The X-ray 
response (in e~) is given by the product of the mea- 
sured 55 Fe peak ADU intensity (measured in Sec- 
tion [2]) with the "true" conversion gain (measured 
in Section [3]). We must first divide the gain by the 
detector internal QE to account for charge that may 
have recombined before being collected in the pix- 



els' depletion regions. Figure [6] plots the measured 
QE of each detector. 

The relevant QE measurement from Figure [S] de- 
pends on the wavelength of the photons with similar 
penetration depths as 55 Fe X-rays. The X-ray pen- 
etration in HgCdTe follows an exponentially falling 
probability distribution. The attenuation length is 
defined as the depth, x, into the material, measured 
along the surface normal, where the intensity of X- 
rays falls to 1/e of its value at the surface. Our 
model, 



2.025 - 9.472 x 10 _ V + 1.8 x 10 



C 7' 



(17) 



is based on data from lHenke et all (|1993f ) and sub- 
sequent modeling 10 , where x is given in microns and 
e 7 is the X-ray energy in eV. The 55 Fe X-ray pene- 
tration depth is towards the very back of these de- 
tectors an d has an associated cutoff wa velength of 
~1.68 (im (jSchacham fc Finkmanlll985h . 

The QE of a detector, however, typically combines 
the effects of several external factors, such as the 
number of photons that pass through the detector 
and the number of photons reflected by the detec- 
tor's anti-reflective (AR) coating, and the internal 
QE, which includes only losses in the detector crys- 
tal. Since our tests measure the detector response to 
individual X-rays, which are either absorbed or not, 
we are concerned only with the internal QE. For the 
plot in Figure O we subtract the AR coating contri- 
bution with a model AR curve. We also assume that 
all photons are absorbed as long as they are short- 
ward of the roll-off at long wavelengths. Therefore, 
instead of measuring the QE at precisely the as- 
sociated cutoff wavelength calculated above (^1.68 
^m), we measure the QE just before the sharp roll- 
off, as indicated by the vertical line. We discuss 
the systematic error associated with this method in 
Section [431 

Table [2] lists the resulting pair-production energy 
for each detector. Figure[7]plots the pair-production 
energy versus cutoff wavelength and bandgap en- 
ergy. FPAs 129 and 153 stand out from the rest. 
Interestingly, these 2 detectors exhibit traits and 
data sets uncharacteristic of the typical FPAs that 
we tested. 

For example, FPA 129 has a large error bar trace- 
able to the unique data available for this FPA. Un- 
like the other FPAs, the conversion gain data for 
FPA 129 consist of only 5 data point s, all of which 
have signal levels >5000 ADU. Figure 5(b) indicates 
that the slope of the conversion gain only begins 
to change at signal levels <5000 ADU. Considering 
the 55 Fe intensity is ~800 ADU, we approximate 
the lower signal gain from an analysis of the photon 
transfer of warm pixels in dark ramps, which has a 
large associated error. The resulting extrapolation 
of the FPA 129 data set is significantly less con- 
strained and, subsequently, has a much larger sys- 
tematic error than the other FPAs. Given the lim- 
ited data set, we do not believe that the error associ- 
ated with FPA 129 represents the typical systematic 

10 http: / /henke.lbl. gov / opticaLconstants /atten2 .html 



error associated with our procedure and therefore 
exclude this data set from our final average. 

The data point for FPA 153 in Figure [7] is clearly 
isolated from the other FPAs. FPA 153 has a sig- 
nificantly lower QE than the other FPAs, especially 
at the cutoff wavelength. In fact, we find that the 
QE for this FPA tends to swing ±20% between dif- 
ferent data sets. Furthermore, FPA 153 shows the 
most rapid decline in conversion gain at low signal 
levels. We do not completely understand the odd 
behavior FPA 153 exhibits, but we recognize that 
this is not standard behavior for these FPAs. We 
therefore do not include the measurement in our fi- 
nal average. 

Averaging over the remaining FPAs, the pair- 
production energy of 1.7 /mi HgCdTe is 2.61 ± 0.16 
eV, which corresponds to 2273 ± 137 electrons pro- 
duced by each 55 Fe X-ray. For comparison purposes, 
Figure [7] plots the pair-production energy of sev- 
eral other common semiconductor m aterials, as well 
as the semi-empirical trend given by Ali g fc Blooml 
(fT975l) . The average pair-production energy for 1.7 
/im HgCdTe falls within a few per cent of the value 
predicted bv lAlig fc Blooml (|1975h . which is within 
the given error. We did not necessarily predict this 
agreement given that the trend-line is only semi- 
empirical. In fact, the measured pair-production 
energy for germanium does not fall e ven close to the 
trend-line. Nonetheless, Figure 1 of I Alig fc Blooml 
(|1975f ) reveals that a majority of materials do fall 
reasonably close to the line. 

The overall consistency of the 1.7 /im results indi- 
cates that the 55 Fe results may be an easy, straight- 
forward sanity check for a given FPA. For exam- 
ple, the inconsistent results for FPAs 129 and 153 
indicated a bad data set and a misbehaving FPA, 
respectively. In other words, if the calculated pair- 
production energy for a given FPA falls far from the 
predicted value, the testers may wish to review the 
other FPA characteristics. 

4.2. Statistical Errors 

We identify four possible statistical errors in our 
measurements: the (1) detector response to the X- 
rays in ADU, (2) the QE, (3) the conversion gain, 
and (4) the coupling coefficient, a. Given the large 
number of samples in all cases, the associated statis- 
tical errors are negligible. Specifically, we determine 
the detector response from the 55 Fe peak in Figure 
[2J In this case, the large number of pixels greatly 
reduces the error of the mean to < 1%. Similarly, 
the large number of pixels reduces the error bars of 
the measured QE and conversion gain to < 1%. The 
curves fit to the gain data subsequently have negli- 
gible error bars. Finally, Equation [TO] indicates that 
a coupling coefficient of ~1% requires > 10 6 sam- 
ples to obtain a signal-to-noise > 10. In terms of 
our autocorrelation matrices, individual exposures 
> 10 6 pixel samples (1024x1024 array) easily sat- 
isfy this requirement. 

4.3. Systematic Errors 

The systematic errors tend to be larger than 
the statistical errors described above. Specifically, 



we identify three possible systematic errors in our 
methods: measurements of the (1) detector re- 
sponse to the X-rays in ADU, (2) the QE, and 
(3) the conversion gain. We measure the detector 
respo nse peak by fitting the histogram in Figure 
2(b) with a combined gaussian and polynomial us- 
mglDL's GAUSSIAN routine. The peak of the dis- 
tribution varies, however, depending on the signal 
thresh old ch osen to indicate a well-centered hit in 
Figure 2(a) We set the threshold at the "kink" 
where the low-intensity plateau and gaussian peak 
meet. This point is rather arbitrary and may not 
be optimal. Varying the threshold by a reason- 
able amount (e.g. <10 ADU) i n cith er direction 



shifts the gaussian peak in Figure [2(b) by a similar 
amount. This shift, which typically results in only 
~1% error, presents only a minor contribution to 
the overall systematic error. 

Measuring the QE introduces an additional sys- 
tematic error. As discussed in Section 14. 1[ the 
55 Fe X-rays have an associated cut-off wavelength of 
~1.68 fim. At these long wavelengths, however, the 
QE is a combination of both internal losses and un- 
absorbed photons. Concerned only about electrons 
generated by absorbed X-rays, we desire a QE that 
measures internal losses only. We assume that all 
photons are absorbed short-ward of the QE roll-off 
at long wavelengths in Figure [5] Instead of measur- 
ing the QE at 1.68 /im, we estimate a shorter wave- 
length at which the QE represents internal losses 
only, as indicated by the vertical line in Figure O 
We assign an error of 0.05 /im to this estimated 
wavelength as it is not entirely clear where the roll- 
off begins. This error corresponds to a systematic 
error for the QE, which depends on the degree of 
slope in the QE curve. For all FPAs, the resulting 
QE systematic error is <5%. 

The measurement of the conversion gain intro- 
duces the largest systematic error to our results. 
Figure [5] shows the measured gain as a function of 
signal level. The desired gain for our purposes must 
be measured at signal levels comparable to the X- 
ray intensities, which are typically ^800 ADU. Our 
archival gain data, however, is limited to higher sig- 
nal levels than necessary. As discussed in Section 
13.41 we therefore choose to fit a curve to the high 
signal gain data and extrapolate to lower signal lev- 
els. Gi ven th e c ontra sting nature of the trends in 
Figures 5(a) and |5(b) we find the likely conversion 
gain at low signal levels is bound by the extremes of 
both an exponential and polynomial fit. We there- 
fore take the conversion gain to be an average of 
these two fits. The associated systematic error is 
given by the difference between the average and the 
values of the two curves at the desired intensity (in 
this case, ~800 ADU). We find a systematic error 
of <10% in all cases. 

The systematic errors are listed in Table [2] and 
included on the plot in Figure While these sys- 
tematic errors contribute to no more than ^11% in 
quadrature, we believe these errors can be signif- 
icantly reduced with more appropriate data. For 
example, we measure the gain by extrapolating our 



data sets down to low signal levels. Future tests 
should be sure to measure the gain directly at these 
signal levels. Such foresight would reduce the sys- 
tematic errors by at least 50%, if not more. 

5. SUMMARY 

In this paper, we have calibrated the 55 Fe X-ray 
response of Teledyne 1.7 /im HgCdTe flight grade 
detectors for the HST WFC3. Our conversion of 
the X-ray intensities from ADU into electrons im- 
plements a technique that restores the "true" gain 
of the detectors using the classical propagation of 
errors. The advantage to this technique is that it 
allows for a direct measure of the covariance ma- 
trix from the data. Monte Carlo simulations have 
validated this technique as being most effective at 
restoring variance that is lost due to charge coupling 
effects. 

The sample of eight 1.7 /im H1R detectors 
yields 2273 ± 137 electrons per 55 Fe X-ray, which 
corresponds to a pair- production energy for 1.7 
/im HgCdTe of 2.61 ± 0.16 eV. This result includes 
both statistical and systematic errors. Better de- 
signed tests can significantly reduce these errors in 
the future. Nonetheless, these results now make the 
55 Fe gain measurement technique of potential inter- 
est for near-infrared arrays, as it has been in the 
CCD community for many years. 

The 55 Fe X-ray technique provides a valid mea- 
surement of the conversion gain at signal levels com- 
parable to the 55 Fe X-ray intensity. This low signal 
level (~800 ADU) is consistent with the signal lev- 
els at which the gain is typically reported for WFC3 
detectors. This technique provides a straightfor- 
ward alternative to the photon transfer technique, 
which can overestimate the conversion gain within 
HgCdTe detectors by up to ~20% (and Si-PIN ar- 
rays by more than a factor of 2) when not account- 
ing for charge coupling effects, such as IPC. The 
known detector response to 55 Fe also will aid future 
detector design and performance studies, such as 
measuring drifts in the system gain or the internal 
QE. 

Although these measurements were made using 
Teledyne HAWAII detectors, we believe that they 
are likely to be applicable to HgCdTe detectors 
from other vendors, so long as differences in cut- 
off wavelength are correctly accounted for. Users 
of Teledyne's older substrate-on should not use this 
metho d as the substrate absorbs the X-rays. Al- 
ready, iFinger et al.l (|2008h have performed similar 
measurements of the 55 Fe X-ray response for 2.5 
/im HgCdTe H2RG detectors. They report results 
consistent with the overall trend in Figure [3 but 
these results fall significantly above the predicted 
line compared to the 1.7 /im results. The authors 
use different calibration techniques and do not re- 
port a detailed review of their systematic errors. In 
the future, we plan to implement our method to 
check the 2.5 /im results, as well as test HgCdTe 
with other cutoff wavelengths (e.g. 5 /im). 

This research was supported by NASA as part of 
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APPENDIX 

A. HOW TO PROCEED WITHOUT ASSUMING THE FORM OF THE COVARIANCE 

MATRIX AS A PRIOR 

A.l. Correlated Noise & Partially Symmetric PSF 

In most HST detector arrays, the detector PSF is symmetric in reflection about the x and y axes to within 
a few a few percent. For the most accurate measurements of conversion gain, this asymmetry needs to be 
accounted for when it is detectable. In this case, the detector PSF can be written as follows, 

/ a \ 

PSF= l-2(a+/3) [3 , (Al) 



a 



and the deconvolution kernel is, 

a(3 



PSF = 



(-1 + 3a) (-1 + 3/3)(-l + 3a + 3/3) 

-2 + 3a + 3/3 -2 + 3a + ±=j^ + 3/3 -2 + 3a + 3/3 

- 2 + 3a + ±^ + 3/3 -l+3a 2 (-l+ffl + (4-y+a(4+/3(-ll+3ffl) _ 2 + 3^ + l^M + 3/ j | (A2 ) 

-2 + 3a + 3/3 -2 + 3a+i^ + 3/3 -2 + 3a + 3/3 

Following the same methods as were used in Section [3+] we arrive at an expression for a\ that is accurate 
to second order in a and /3, 

$4 = C*4,4 + a (— Ci t 4 — C^i + 4C4_4 — C?4,7 ~ CV,4) 

+ /3 ( — C 3j 4 — C*4 :3 + 46*4,4 — 64,5 — 6*5,4) 

+ a 2 (Ci.x - 5C M + Ci, 7 - 5C 4 ,i + 16C 4 , 4 - 5C 4 , 7 + C 7 ,i - 5C 7 , 4 + C 7 , 7 ) 
+ /? 2 (C 3 ,3 - 5C 3 ,4 + C 3 , 5 - 5C 4 , 3 + 16C 4 , 4 - 5C 4 , 5 + C 5 , 3 - 5C* 5 , 4 + C 5 ,b) 
+ a/3 (2C ,4 + Ci,a - 6C M + Ci, 8 + 2C 2 , 4 + Os,i - 6C 3 , 4 + C 3 , 7 + 2C 4 , - 6C 4 ,i + 2C 4 , 2 

— 6C*4, 3 + 24C4,4 — 6(^4,5 + 2C4 ; 6 — 6C4.7 + 2C4.g + C*5,l — 6(75,4 + Cs :7 + 2C*6 l4 

+C 7 , 3 - 6CV, 4 + C 7 , 5 + 2C 8 , 4 ) (A3) 

Because the mathematics are complex, we do not show the steps here, other than to note that all computa- 
tions were done using the Mathematica symbolic mathematics program. 

A. 2. Correlated Noise & Completely Symmetric PSF 

For our calibration of 55 Fe X-rays, we believe that the method described in this paper (with extensions 
to remove the need to assume a functional form for the detector PSF as a prior) is likely to have the 
advantage that the covariance matrix can be explicitly computed from the data to serve as a diagnostic of 
other correlations that may be present in the data in addition to IPC. Unfortunately, doing this requires 
specialized detector readout modes that have not been implemented for HST because they are not required 
by the planned science. 

If we do not assume a theoretical form for the covariance matrix, Equation |T2] can be written as follows, 

d\ « 6*4,4 + a (— Ci,4 — C 3j 4 — C 4) i — C*4, 3 + 864,4 — 64,5 — C*4 j7 — 65,4 — 67,4) 

+ a 2 (2C*o,4 + C*i : i + C*i, 3 — 11C*i : 4 + C1.5 + C*i : 7 + 2C*2,4 + C 3 ,i + C 3 , 3 — 11C 3 ,4 + 6*3,5 
+C* 3 ,7 + 2C 4 , - 11C!m + 2C 4,2 - HC 4 , 3 + 56C 4 , 4 - 11C 4 , 5 + 2C 4 , 6 - 116*4,7 + 2C 4 , 8 
+C* 5 ,i + Cr , 3 - 11C M + C* 5 ,5 + C 5 , 7 + 2C 6 , 4 + CV.i + C 7 , 3 - l\C 7A + C 7 , 5 + C* 7 , 7 + 2C* 8 , 4 ) + . . . . (A4) 
The individual C{j can then be measured using the relation, 

^ n—l 

C i,3 = ( Sl ~ ( S 3 ~ ?J)„ ■ ( A5 ) 

k=0 

In Equation lA5[ n refers to the exposure sequence number, which needs to be large for accurate measurement. 
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TABLE 1 
Autocorrelation Matrix for 
FPA 152 



1.000 


0.041 


0.006 


0.005 


0.021 


0.005 


0.001 


0.001 


1.4e-4 


7.5e-4 


2.6e-4 


-8.0e-6 


-3.5e-4 


0.001 


7.1c-4 


-4.0c-4 



TABLE 2 

Summary of Detector Properties 



Serial Number ov ca< j c \ / TT nr ^ i ^ ±- \ Gain Correction pair-production Energy 

(bctorc correction) (Hor/vcrt) (after correction) c e oj 





e rms 


e-/ADU 


% 


e-/ADU 


% 


eV 


FPA 129 


13.74 


3.25 


2.98/1.50 


2.55 


21.6 


2.52 (0.012 ± 0.31) 


FPA 148 


14.59 


2.64 


1.93/1.17 


2.34 


11.6 


2.64 (0.002 ± 0.07) 


FPA 150 


19.80 


2.73 


2.07/1.29 


2.38 


12.7 


2.65 (0.002 ± 0.10) 


FPA 152 


19.04 


2.72 


1.93/1.12 


2.40 


11.9 


2.47 (0.005 ± 0.03) 


FPA 153 


18.02 


2.67 


2.02/1.11 


2.35 


11.9 


1.98 (0.005 ± 0.09) 


FPA 154 


18.38 


2.67 


1.97/1.25 


2.23 


16.6 


2.60 (0.003 ± 0.05) 


FPA 160 


15.68 


2.67 


2.04/1.28 


2.30 


14.0 


2.64 (0.002 ± 0.04) 


FPA 165 


13.37 


2.62 


2.22/1.46 


2.18 


17.0 


2.68 (0.007 ± 0.06) 




(a) Symmetric PSF 





a 






1-2(a+|3) 






a 





(b) Partially Symmetric PSF 



Fig. 1. — Illustration of charge coupling mechanism for both the (a) symmetric and (b) partially-symmetric cases. In both 
cases, photocurrent entering a central pixel may be displaced into neighboring pixels by coupling mechanisms such as interpixcl 
capacitance (IPC). The fraction of charge coupled into the neighboring pixels is given by the coupling coefficients, a and j3. 
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(b) Nearest Neighbor Summation 

Fig , 2. — Histograms of the observed 55 Fe X-ray response, in AD U. for a typical 1.7 /im HgCdTe detector (FPA 152). Figure 
|2(a)| shows only the individual pixel response, while Figure |2(b)| includes the summation of the nearest neighbor intensities. 
The Kq intensity peak is sharper and better defined by the individual pixels, but the histogram ignores charge coupling 
effects. Adding the nearest neighbor pixel intensities compensates for charge coupling and thereby provides a more accurate 
measurement of the true X-ray response. 
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Fractional Coupling (%) 

Fig. 3. — Histogram of fractional crosstalk observed in neighboring pixels when compared to the total X-ray intensity of 
a typical detector (FPA 152). The nearest neighbors show a significant degree (1-2%) of fractional crosstalk with the central 
pixel. In contrast, a distribution of random pixels shows no coupling the X-rays. 
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Fig. 4. — Monte Carlo validation of Equation 1 131 Of the models that were tested, the propagation of errors technique 
best recovers the "true" variance. The four models that were tested include: no correlated noise correction ("No IPC corr."), 
multiplication by corrective coefficient ("1 + 8a Coefficient"), implementation of equation 34 from Moore ct al. (2006) ("Auto- 
correlation"), and application of the correlated noise correction using propagation of errors ("Prop, of Errors"). 
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(a) FPA 152 
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(b) FPA 153 



Fig. 5. — Gain (e~/ADU) versus signal intensity (ADU) for FPAs 152 and 153. The corrected gain has had the "true" 
variance restored via the propagation of errors. As predicted by the Monte Carlo simulations, the uncorrected gain is typically 
overestimated by ~ 13%. The relevant gain is measured at signal levels consistent with the 55 Fe X-ray intensity, which are 
typically ~800-1000 ADU, as indicated by the vertical lines. We measure this number by extrapolating a best fitting line down 
to lower signal levels. While both the polynomial (dotted line) and exponential (solid line) fits yield similar results in the case 
of FPA 152 (a), the polynomial fit appears to be the most appropriate fit in cases where the gain rapidly declines at lower signal 
levels, as in the case of FPA 153 (b). (FPA 153 is the worst case within our data set.) We therefore prefer the polynomial fit 
for our extrapolation in all cases. 
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Fig. 6. — Measured quantum efficiency (QE) of each detector, as well as the reflective coating transmission model, as a 
function of wavelength. The vertical line denotes the photon wavelength corresponding to the equivalent penetration depth 
of the 55 Fe X-rays. Like these long wavelength photons, the X-rays typically penetrate deep into the detector. We take the 
equivalent QE at longer wavelengths just before the sharp roll-off, at which the electron capture efficiency no longer dominates 
the QE measurements. 




Fig. 7. — Measured pair creation energies, e, of the HST WFC3 1.7 ;im HgCdTe detectors plotted against the semiconductor 
bandgap energy, E g , and the cutoff wave length correspond ing to E g , A g = hcEg~ . Also plotted in this figure are several other 
semiconductor materials published in Alig & Bloom (1975)'s figure 1. The solid line is the semi-empirical relationship between 
e and E g , derived in I Alig fc Blooml (197a l. 



